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On Ignoring the Singularity in Numerical Quadrature. 

B/ 

R.K. Miller 

Introduction 

Davis and Rabinswitz [l] recently studied the question 
of "ignoring the singularity " in numerical quadrature. 
That is, if f(t) becomes singular at a point £ where 
a £ £ <^ b then one defines f(£) = 0 (or any other 
finite value) and then approximates the integral. 

I = / b f(t)dt 
a 

by a usual numerical quadrature rule. They show that 
this procedure is not valid in general. However if 
£ = a (or some other rational point of [a,b]), then 
compound quadrature rules do approximate I when f is 
monotone near £. Certain of the positive results in 
M where generalizied by Rabinswitz 0] . Gautschi 
[3] has applied a result in [2] to two other quadratures 
of interpolatory type. 

The purpose of this paper is to generalize some of the 
convergence theorems in [i] and [2]- We shall replace 
the assumption of monotonicity of f(t) near t = £ by the 
more general condition that f(t) can be dominated by a 
monotone, integrable function. We shall also establish 
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some theorems on error bounds and convergence rates which 
are similar to those obtained in !>]. 

2. General Quadratures 
Let M be the set 

M = {fEC(0 > T]rU 1 (0,T) : f is nonnegative and non- 

increasing on 0 < t 5 T} . 

Define M to be the set of all functions f £ C(0,Tl 
d J 

such that f can be majorizied ty a function in M, 

M d = {f E C(0,T]:3 FEM with | f (t ) |sF(t ) on 0 < 
t 2S T } . 

For any function f E M assign the arbitrary value 

d 

f(t) = 0 when t = 0. 

Given any numerical quadrature rule 

Q(f) = Ej =0 W(J) f ( T ( j ) ) , f E C[0,T] 

let Q(f,S) be the modified quadrature rule obtained from 
Q by redefining f(0) = 0, 

Q(f,S) = £j= 0 


W ( j )f (T( j ) ) , f (0) := 0 
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Then Q(f,s) ignores possible singularities at t = 0 
and is well defined for all functions f £ M^. In general 
Q(f,S) 7 1 Q (f) for all functions f£C[o,T3* With these 
preliiainaries we are ready to generalize the lemma 0 ]. 
First consider rules which are open at t = 0 . 

Lemma 1. Consider a sequence of rules 

V f > ■ VJ> f(T n ( J>> 


where 


0 < T ( 0 ) < T (l) < ... <T (m ) < T 
n n n n — 


an d W (j)>0 for all j . Define T q (-1) = 0. 

Suppose there exist positive constants C and A such 

that uniformly for all positive inte ger s n , i_f 

j = 0(l)m and if |T (.1)1 < A then 
n 1 n 1 

(2-1) W n (j) < C (T n (j) - T n (j-1)}. 

(j = 0 (l)m n means j = 0 , 1 ,2 , . . . ,m n . ) 


V 
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Suppose for each function g E C [ 0 , T] one 


(2-2) lim Q n (e) = /g3(t )dt , n->°° . 


Then for any function fEM, 
a 


(2-3) lim Q n (f,S) = /Jf(t)dt, n-°°. 


In particular if 0<B<A, if one defines 


f B (t) = f(t) on B < t < T; = f(B) or 


and i f 


6(t) = sup { | f ( s ) - f B ( s ) | : t < s < 


then the error 


E s ( f , Q n ) = /q f(t)dt - Q n (f,S) 


satisfies the estimate 


(2-1*) | E s ( f ,Q n ) | < |E(f B ,Q n ) | + | f 


f ( B ) }dt | + C 6 ( t )dt . 


has 


0 < t < B, 


T), 


(f(t) - 
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Proof. Write E = E g (f,Q n ) in the form 

E = /q{ f ( -t ) - f B (t)}dt + E(f B> Q n ) + W n (j) 

{ f ( T ( j ) ) - f(T (J)} 
i) n n 

= /®{r(t) - f(B)}4t + E(f B ,Q n ) + e n< 


Then for any n 


£ 

n 


< Z 


m 

n 


j = 0 


W n (j) l f < T n (j)) - f B (T n (j)) l 


< 


Z 


m 

n 


j = 0 


W n (j) 6 (T n (j)) = Q n (6, S) 


Since f£M, , there exists a majorizing function F£M. 
a 

Then for s in the range b<t<s<B one has 


| f ( s ) - f B (s)| < F(s) + F(B) < 2F(t). 


Therefore 


6 (t) < , 2F(t) on 0 < t < B, 6(t) = 0 

on B £t<l , and hence 66L^(0,l). Note also that 6(t) 
is continuous, nonnegative and noni ncr e as ing . This 
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together with (2.1) implies that 

Q n (6 ’ S) = £ T n (j)<B 1 CI{T „ ( J> - 

T n (J-l)}6(T n (j) < C I /*»{”!, 4 

( t ) dt = C /® 6 (t)dt. 


This completes the proof of (2 • ^ . 

Line (2.3) follows immediately from (2.U) and the 
estimate 6 ( t ) £ 2 F(t). Indeed by first chocsing B 
small and then choosing n large one can make the right 
hand side of (2.b) as small as desired. Q.E.D. 

Almost the same result is true for quadratures which are 
closed at t = 0 . 


Lemma 2 . Suppose a sequence of quadrature rules Q 


satisfies the two conditions 


0 = T (0) < T (l) < ... < T (m ) < T 
n n n n — 


w n (j )>0. 


Suppose there exist positive constants C and A such that 
if | ( J ) | < A and if j = l(l)m n then (2.1) is true 
uniformly in n. I_f (2.2) is also true then (2.3) follow s . 
In particular if f, E C[0,T] and 66M a are the functions 

* 5 £ 

defined in Lemma 1 and i f 0 < B < A, then 
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(2.5) |E s (f,Q n )| < | Q n ) | + |/*{f(t) - 

f ( K ) } d t | + C 6 ( t ) dt + W (0) |f(B)|. 

Proof . The proof is the same as that of Lemma 1 except 

for the estimates of e . In this case 

n 

e n * Vo W n ( J> - V T n (J)1 - 


and 


|e n l < Q n (5,S) + W n (0) | f (B) | 

< C/®6(t)dt + W n (0) | f (B) | 


vhere we define 6(0) = 0. Our hypotheses easily imply 

that W (0) -*■ 0 as n -*■ °°. Therefore Lemma 2 is proved, 
n r 

Combining the two results we have proved: 

Theorem 1. Consider a sequence of numerical quadrature 
rules Q where 


0 < T (0) < T (1) <...< T (m ) < T and W (j)>0. 
— n n n n — n 


Suppose (2.1) is true as in Lemma 1 (when T ( 0 ) > 0 ) 
or Le m ma 2 ( T ( 0 ) = 0 ) . I_f (2.2) is also true , then for 
any f E 



e 


Q ( f ) *» s' 1 ' f (t ) it as n -► 
n u 

Indeed if 0 < B < A end if f n € C [0,T] and 6 G M, 
£ 

are the fun ctions define d i_n Lemm a 1 t h c n 

(2-6) |E s (f,Q n )| < |E(f B ,Q n )| + |/J(f(t) - 

f ( B ) } dt | + C/® fi(t)dt + {1 - s C n T n (0)} 
W n (0) I f ( B ) I . 

One can also generalize the Corollary in [2, p. 196 ] in 
the obvious way. 


3. Compound Rules. 

Consider a quadrature rule R defined on the interval 
0 < t < 1 : 


(R) R(t) = 


Vo V'V 


where J > 0 and 


(3-1) 0 < t„ < t, < ... < tj _ 1 , Wj > 0 , ij = 0 W = i 


- 0 1 


(If t^ > 0 then define t 


= 0.) For any integer n 1 
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ana any interval 0 <_ t <_ T one can then define a 
compound rule 

(n X R) R n (f) = = f(t J H + Kk)} 

where H = T/n. hot C>0 be any constant satisfying 
(3.2) Wj < (tj - t j _ 1 ) C (J = l(l)J) 

and either 

(3.3a) W Q < (t Q + 1 - tj) C (if t Q > 0 or tj < l) 

or 

(3.3b) (W Q + Wj) X (tj - t J _ 1 ) C (if t Q =0 and tj=l). 

Theorem 2. I_f (R) satisfies (3.1) then for any f E 

lim R (f,S) = /^f ( t ) dt . 
n u 

n->°° 

Proof . The definition (3.2-3) of C implies that (2.1) 
is true with A=T. Since R integrates constants and 
n -*°o then (2.2) is also trival. Therefore Theorem 2 is 


a corollary of Theorem 1. 


Q . E . D . 
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The error estimate ( 2 . 6 ) is rather pessimistic for 
compound rules . Therefore we shall derive another 
estimate which is more suitable for many purposes. 

Let K>0 be the smallest constant which satisfies 

(3.I.) Vj < (♦, ■ t,^) K 

for j = 1(1)J if t Q =0 or for j = 0(l)j if t >0 . 

Theorem 3. Suppose (3.1) and (3.^) are satisfied . Let 

H = T/n for any function f GM, define 
a 

f u (t) = f ( H ) if 0 < t < H; f(t) if H < t < T. 

n 

I f FGM is any majorizing function for f then 

t H 

(3.5) |E s (f,R n )| < |E(f H ,R n )| + /"F(t)dt + K F ( t )dt 
or 


l^ f » R m )| 1 l E(f H’ R n ) l + (1 + K) ; 0 F(t)dt * 

Proof . Since (R) integrates constants, then the error 
may be written in the form 

E s (f ’ R n> “ E(f H> R „> + E 0 
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where 

(3.6) E q = f Q f ( t ) dt _ HVjf ( t ,H) 

Since |f(t)| < F(t) cn 0 < t <_ T, then 

I/® f (t)dt ] < r(t )dt . 

Let a = 0 if t Q > 0 and a = 1 if t Q =0. 
and ( 3 . ^ ) imply 

| £j = 0 H Vj f ( t , H ) | < l J a Hw^ F(t,H) 

< K E J H(t - t ) F(t,H) 

— a j J -1 J 

T t • ^ -ft 

< Z J / J F(t)dt = K / n 

” ° *j-i 0 

This proves (3-5) and the theorem. 

If one knows that fEC 1 (0,T], then the 
may be estimated using Peano's theorem. 
That is 


B'W = 'l P n <S) f i (s)dt 


Then f(0) = 0 


F ( t ) dt . 


term E(f u ,R ) 
n n 


where P q is the appropriate Pear.o kernal . 


Since 
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P (s + H) = P(s)onO«-s< T-H, then I? (s)| need be 
be estimated only on the interval 0 < s < H. Therefore 
the following result is an immediate corollary of Theorem 3. 

Corollary 1. Assume the hypotheses of Theorem 3. If. 
f G C 1 (0,T] t hen 

(3.7) | E(f »R n ) | < | |P n | | /£|f'(t)|dt + /J F(t)dt + 

t H 

k ; o m F(t)dt 

where ll p n ll = SU P {|P n (s)| = 0 < s < H). 

Corollary 1 is useful in estimating convergence rates 
in certain c as es • Foil owing [l|] we shall say that a 
function f E C^(0,T] is weakly singular at t = 0 if the 
function 

a(t,f) = | f (t ) | + I’l |f ' (s) |ds 
is in ( 0 , T ) . 

Corollary 2. Suppose the hypotheses of Theorem 3 are true . 
If f is weakly singular ( a_t t = 0 ) t hen 

(3.8) E ( f , R ) = 0(/n a (t ,f )dt ) as H+0 . 

s n u — 
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Proof . It follows immediately from the definition of 
weakly singular functions that f E and tx(t,f) E M 
is a majorizing function. Thus (3.7) implies 

l E s ( f > R n H 1 I I p n I I /j5|f'(t)dt + a(t,f)dt + 

K/|Ja(t ,f )dt 

< | |P n | | a(H,f) + (1 + K) a(t,f)dt. 

Using the estimate | | | | 2H (see fpr example [ k , 

section II ] ) and the monotonicity of o. one has 

| |P n I I «(H,f) < 2H o(H,f)< 2/“a(t,f)dt. 

Therefore for any n (H = T/n) one has 

| E( f ,R n ) | < ( 3 + K) /Jj a(t ,f)dt. 

Q.E.D. 

For example if f(t) = t (0<p<l) then (3.8) predicts 

that E (f,R ) is at least of order h* If f(t) = t ^ 

s n 

sin (t 3.) where 0 <p,q <1 and p + q <1 then E (f,R ) is 

s n 

at least of order h^ ^ If f(t) = t ^ sin (t where 

0 < P < 1, q>0 and P + q > 1 then our theory predicts 
convergence but gives no order estimate. 
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4. Numerical Example. 

The data in [l] will he used to illustrate the theory 
given above. For the midpoint rule M(f) = f(|) one has 
H = h = T/n. Since -P^Cs) = s if 0 < s < h/2; 2-h if 
h/2 < s < h, then | |p | | = h/2 and K = 2. Therefore (3.2) 
has the form 

|E s (f,M n )| < (h/2) /^ | f * ( t ) | dt + /q F ( t ) dt + 2/£ /2 F(t)dt. 
Table 1 contains data for the case 
(4.1) /J t" l/2 dt = 2. 

The fourth column is the theoretical error computed using 
( 3 . T ) • This error bound is seen to be pessimistic by a 
factor of 7 to 8. Corollary 2 suggests that the error may 
be of the approximate form 


(4.2) E (f,M ) = C/h 
s * n 


(h=T/n) 


for some constant C > 0. The ratios E (f,M ) / E (f.M ) 

s n s' n+1 

are given in column five. The theoretical ration computed 
using (4.2) is /2 (column six). It can be seen that (4.2) 
is approximately true with C = .6l. 


(insert table 1 near here) 
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Table 2 contains similar data for ( l* . 1 ) using the 

trapezoid rule ( T) , Simpson's rule (S ) and the Gaussian 

two point rule (G^). The theoretical errors are good 

for the trapezoid rule and progressive worse for Simpson 

and Gauss two point. In all cases [k . 2 ) is approximately 

true, E ( f , T ) = 1.5/h, E (f,S ) = . 89 /h and E (f,nXG 0 ) 

’ s n s n s 2 

= .35/h. 

(insert table 2 near here). 

One could also analysis the data in [l] for ( . 1 ) using 
the Gauss U8 - point rule. In this case E g (f,n X G^g) 

= .l8/h (approximately). It would be very difficult to 
estimate | | P | | accurately in this case. Thus (3.7) is 
essentially useless. 

Data for the example 

1 1 1 

t~? sin (t“ ? ) dt = 1.0813U 

is given in table 3 for Simpson's rule and for the Gauss 
1*8 - point rule. For this case our . theory predicts conver- 
gence but yeilds no useful information on errors or conver- 
gence rates. For Simpson's rule the method appears to be 
convergent. For G^g the method starts to converge nicely 
but blows up at the fourth step. 


TABLE 1 



M n (f 0> 

Error 

Th .Error 

Ratio 

Th . Rat i 0 

5 

1.8931 

.1069 

.6763 

1 . 1* 1 1* 

1 . l*lU 

6 

1.921*1* 

.0756 

.1*815 

1.1*13 

l.l*ll* 

7 

1.9^65 

• 0535 

.1*321 

1 .1*15 

If 

8 

1.9622 

.0378 

.21*27 

1.1*16 

If 

9 

1.9733 

.0267 

.1720 

1 .1*13 

II 

,10 

1.9811 

.0181* 

.1218 

1 .1*10 

ft 

,11 

1.9866 

.0131+ 

.0861* 





TABLE 2 


2 5 X 

T 

Approx . 
1.7418 

Er ro 
.2582 

2 6 X 

T 

1.8174 

.1826 

2 7 X 

T 

1.8709 

.1291 

X 

co 

CVJ 

T 

1.9087 

.0913 

2 9 X 

T 

1.9355 

.0645 

X 
o 
1 — 1 
OJ 

T 

1.9544 

. 0 u 56 

2 5 X 

S 

1.8427 

.1573 

2 6 X 

s 

1.8887 

.1113 

2 7 X 

s 

1.9213 

.0787 

2 8 X 

s 

1 . 9444 

• 0556 

2 9 X 

s 

1.9606 

.0394 

2 10 X 

s 

1.9721 

.0297 

2 X 

G 2 

1.7528 

.2472 

It X 

G 2 

1.8252 

.1748 

6 X 

°2 

1.8573 

.1427 

8 X 

G 2 

1.8764 

.1236 

10 X 

°2 

1.8894 

.1106 

12 X 

°2 

1.8991 

.1009 

lU X 

°2 

1.9066 

.0934 

16 x 

°2 

1.9126 

.0874 


Th . Error 

Ratio 

Th . Rat 

.6031 

l.4l4 

1 . 4l4 

.4294 

l.4l4 

1? 

• 3055 

l .4i4 

ft 

.2168 

1.415 

It 

• 1537 

1 . 4l4 

ft 

. 1089 

1.1792 

1.413 

l.4l4 

.8735 

1.4i4 

tf 

.6198 

1.415 

ff 

.4393 

1.411 

ft 

-.3113 

1.412 

ft 

.2203 

7.6572 

1 . 4l4 

l .4i4 

5.4012 

1.225 

1.225 

4.4213 

1.155 

1.155 

3.8341 

1.118 

1.118 

3.4328 

1.096 

1.095 

3.1360 

1.080 

1.080 

2.9050 

1.069 

1.069 

1.7640 



TABLE 3 


2 5 X S 

Approx 

1 .1234 

Er r or 

-.0421 

Ratio 

.241 

2 6 X S 

.9116 

.1697 

1.562 

2 7 X S 

• 9727 

.1086 

1.017 

2 8 X S 

. 97^5 

.1068 

1.745 

2 9 X S 

1.0201 

.0612 

2.696 

2 10 X S 

1.0586 

.0227 


1 X G 48 

. 9UU 9 

.1346 

1.534 

2 X G U 8 

. 9 9 2 

.0889 

2.285 

3 X G U 8 

1.0U02 

.0389 

.2571 

14 x g 48 

.9300 

.1513 
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